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Abstract 


A  novel  theory  for  the  propagation  of  a  spherical  detonation 
wave  is  presented.  The  model  is  that  of  a  point  spherical  blast  in  a 
detonating  medium.  The  reaction  rates  are  assumed  to  be  finite  but 
the  coupling  between  the  shock  motion  and  the  chemical  reactions  is  mo¬ 
delled  by  a  global  function  which  leads  to  an  effective  chemical  energy 
release  at  the  detonation  front  which  depends  on  the  local  shock  strength 
and  the  shock  radius.  The  important  result  obtained  is  the  demonstration 
of  the  existence  of  finite  steady-state  velocities  which  are  below  the 
theoretical  planar  Chapman-Jouguet  value.  The  magnitude  of  these  final 
steady-state  velocities  depends  on  the  value  of  the  ignition  energy 
used  to  initiate  the  detonation  wave  itself.  A  quantitative  comparison 
with  the  recent  velocity  measurements  made  by  Brossard  in  stoichiometric 
and  equi-molar  mixtures  at  .1  ^ 

ment,  lending  support  to  the  present  proposed  theory.  Using  the  experi¬ 
mental  data  of  White  for  the  induction  time,  the  present  theory  can  be 
used  to  predict  the  critical  energy  for  direct  initiation  as  well  as  tne 
transverse  wave  spacing,  and  both  are  found  to  be  in  good  agreement  with 
the  experimental  data. 
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.28  bar  showed  good  agree- 
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1  •  Introduction 

From  the, theoretical  point  of  view,  the  existence  of  Chap- 
man-Jouguet  (C-J)  spherical  detonation  waves  depends  on  the  existence 
of  a  physically  acceptable  solution  to  the  unsteady  gasdynamic  equations 
that  satisfies  the  C-J  boundary  conditions  at  the  detonation  front.  The 
early  studies  by  Jouguet  Taylor  ^  and  Zeldovich  ^  (J-T-Z)  were 

based  on  the  idealized  model  of  a  discontinuous  C-J  front  originating 
at  the  center  f  -  O  at  time  =  O  and  subsequently  expanding  outwards 
at  a  constant  velocity.  The  similarity  solution  they  obtained  possesses 
a  singularity  in  the  form  of  an  infinite  expansion  gradient  immediately 
behind  the  detonation  front.  The  presence  of  this  singularity  throws 
some  doubts  as  to  the  validity  of  the  J-T-Z  solution  on  a  physical  basis, 
hence  the  existence  of  spherical  C-J  detonations  i.s  not  established  by 
the  J-T-Z  model. 

Recognizing  the  experimental  fact  that  the  direct  initiation 

of  a  spherical  detonation  wave  invariably  demands  the  explosive  release 

of  a  finite  quantity  of  energy,  the  reacting  blast  wave  model  (R-B-W) 

(4) 

was  subsequently  proposed  and  analysed  by  Lee  .  The  detonation  front 
is  still  treated  as  a  gasdynamic  discontinuity  but  with  the  inclusion 
of  a  finite  initiation  energy  Sc  no  "teadv  wave  solution  is  possible. 

■  The  detonation  front  is  initially  highly  overdriven  and  propagates 
essentially  as  a  strong  point  blast  wave  (i.e.,  Kg  t  '  ).  As 

the  wave  expands,  the  effect  of  chemical  energy  release  at  the  front 
progressively  influences  the  wave  motion  until  a  C-J  wave  is  approached 
asymptotically  as  /?$  :>  oo  .  A  later  study  of  the  R-B-W  model  by 


Korobeinikov  '  using  a  perturbation  method  demonstrated  a  similar 
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behaviour.  Hence  according  to  the  R-B-W  model,  no  steady  state 
spherical  C-J  detonation  waves  are  possible. 

In  a  recent  study  by  Levin  and  Chernyi  ^  on  the  asymp¬ 
totic  motion  of  a  decaying  c ^erdriven  spherical  detonation,  it  is 
found  that  unlike  the  decay  of  a  non-reacting  blast,  the  final  steady 
state  C-J  condition  can  be  achieved  at  a  finite  radius.  A  more  care¬ 
ful  subsequent  analysis  by  the  Authors  ^  indicates  that  although 
Levin  and  Chernyi's  solution  is  possible,  there  are  no  convincing 
arguments  on  a  physical  basis  in  the  assumption  of  the  particular 
asymptotic  behaviour  of  the  flow  structure  which  is  necessary  to  ar¬ 
rive  at  their  result.  Also,  Levin  and  Chernyi's  analysis  assumes  the 
J-T-Z  similarity  solution  when  the  C-J  condition  is  reached.  From 
the  solution  for  the  R-B-W  model  there  seems  to  he  no  indication  that 
it  would  approach  the  J-T-Z  similarity  solution  as  the  decaying  over¬ 
driven  wave  becomes  a  C-J  wave.  An  attempt  to  verify  Levin  and  Chern¬ 
yi's  asymptotic  solution  by  carrying  out  a  numerical  characteristic 
analysis  of  the  R-B-W  model  to  the  asymptotic  regime  is  inconclusive 

/ON 

as  instability  of  the  type  obtained  by  Fickett  and  Wood  and  Streh- 

(9) 

low  and  Hartung  occurs  when  the  C-J  conditions  are  approached. 

Hence  no  definite  conclusions  can  be  drawn  at  this  stage.  However, 
even  if  the  ■  ,vin  and  Chernyi  solution  does  provide  the  asymptotic  mo¬ 
tion  of  the  detonation  wave  in  the  R-B-W  model,  the  physically  quest¬ 
ionable  J-T-Z  similarity  solution  is  still  retained. 

In  the  J-T-Z  and  R-B-W  models,  the  chemical  reaction  rates 
have  been  assumed  to  be  sufficiently  fast  so  that  the  detonation  front 
can  be  considered  as  a  gasdynamic  discontinuity.  In  the  recent  studies 
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(10)  (1>\ 
by  the  Authors  '  and  Lundstrom  and  Oppenheim  it  is  demonstra¬ 

ted  that  the  rapid  gasdynamic  expansion  immediately  behind  a  spherical 
blast  wave  can  result  in  a  significant  increase  in  the  reaction  zone 
thickness.  Hence  it  is  quite  evident  that  the  discontinuity  assumption 
for  the  detonation  front  would  have  to  be  abandoned  in  a  realistic  theo¬ 
retical  model.  Unfortunately,  a  finite  kinetic  rate  model  where  strong 
coupling  exists  between  the  hydrodynamics  and  the  chemical  reactions 

is  extremely  complex  mathematically.  The  numerical  studies  thus  far  b j 
(12)  (13) 

Rajan  and  Taki  and  Fujiwara  provide  only  a  partial  description 

of  the  problem  in  the  overdriven  regime  of  the  shock  motion.  Instability 
problems  near  the  C-J  regime  of  the  shock  motion  would  be  the  major  dif¬ 
ficulty  in  any  numerical  solution  of  the  exact  model.  It  is  also  doubt¬ 
ful  that  shock  tube  kinetic  data  can  be  used  in  these  exact  models  for 

spherical  detonations  since  the  molecular  energy  transfer  processes  occur 

(14) 

under  highly  non-equilibrium  conditions  .  In  an  attempt  to  reproduce 
qualitatively  the  experimental  observations  for  the  propagation  of  a 
spherical  detonation  wave,  a  simple  phenomenological  model  was  proposed 
by  the  Authors  in  a  recent  paper  ^  \  There  the  coupling  between  the 
chemical  energy  release  in  the  flow  structure  and  the  motion  of  the 
expanding  detonation  front  is  assumed  to  depend  on  the  local  induction 
zone  thickness.  The  dependence  of  the  induction  zone  thickness  on  the 
local  shock  strength  and  the  flow  gradient  is  modelled  globally  by  an 
appropriate  mathematical  functions.  The  results  of  that  preliminary 
study  are  extremely  encouraging  in  that  almost  all  the  essential  features 
of  the  propagation  of  a  spherical  detonation  wave  observed  experimentally 
are  reproduced  on  a  qualitative  basis.  The  present  paper  describes  our 
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subsequent  work  on  this  phenomenological  theory.  Adopting  essentially 
the  same  approach,  the  previous  model  has  been  completely  reformulated 
and  analysed  in  the  present  paper.  The  results  of  the  present  model 
can  now  be  compared  with  experiments  on  a  quantitative  basis. 

2 .  Formulat ion 

As  in  the  R-B-W  model,  we  consider  the  propagation  of  a  point 
spherical  blast  wave  in  a  detonating  gas.  The  blast  is  generated  by 
the  instantaneous  deposition  of  a  finite  quantity  of  energy  at  the 
origin  at  timp  -  O  •  However,  unlike  the  J-T-Z  or  the  R-B-W  model 
the  chemical  kinetic  rates  are  now  considered  to  be  finite  so  that  at 
any  instant  of  time  t  when  the  shock  is  at  a  radius  ,  the 

reaction  front  is  at  (•£)-•  )  lagging  behind  the  shock  by 

a  distance  .  The  present  model  does  not  assume  a  weak  coupling 

between  the  hydrodynamic  and  the  chemical  reactions.  Therefore  the  des¬ 
cription  of  the  shock  motion  and  the  reaction  kinetics  in  the  flow  behind 
the  shock  cannot  be  determined  independently  as  in  the  previous  studies 
of  the  finite  kinetic  rate  blast  model  ^  ^ .  However,  instead  of 

seeking  a  complete  solution  to  the  problem,  which  would  require  the 
simultaneous  solution  of  the  non-linear  hydrodynamic  equations  and  the 
appropriate  kinetic  rate  equations,  we  are  interested  in  the  motion  of 
the  detonation  wave  only.  In  other  words,  we  arc  not  concerned  with 
the  determination  of  the  detailed  hydrodynamic  Clow  structure  and  the 
distributions  of  the  various  chemical  species  in  the  wake  of  the  sphe¬ 
rical  detonation  front,  but  only  in  the  shock  trajectory  a,K* 

the  variation  of  the  strength  of  the  detonation  wave  with  radius.  Within 
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the  context  of  this  modest  aim,  certain  s implifying  assumptions  can 
be  made,  resulting  in  a  cons iderable  reduction  in  mathematical  dif¬ 
ficulties  . 

First  of  all,  we  note  from  the  well-known  blast  wave  theory 
that  the  shock  trajectory  Rstt)  is  rather  insensitive  to  the 
manner  in  which  the  flow  properties  (i.e.,  pressure  j>  ,  density 
and  particle  velocity  (JL  )  are  distributed  behind  the  shock  provided 
that  the  global  conservation  laws  are  satisfied.  By  global  conserva¬ 
tion  laws,  we  mean  the  integral  relationships  for  the  conservation  of 
mass,  i..e., 


(2.1) 


and  the  conservation  of  energy 


where  Q  and  ^Pdenote  the  internal  thermal  energy  and  the  chemical 
energy  per  unit  mass  respectively,  and  the  subscript  "  O  “  represents 
the  initial  undisturbed  condition.  Since  the  flow  is  spherically  sym¬ 
metric,  the  total  momentum  is  always  zero  and  is  conserved  automatically. 
For  blast  waves,  the  most  important  relationship  is  the  energy  integral. 
The  insensitivity  oi  the  shock  trajectory  &<$£&)  to  the  f-ow  struc¬ 
ture  is  noticed  from  the  fact  that  ■'  though  t he  flow  distributions  diifer 
significantly  in  the  various  non-similar  solutions  for  the  blast  wave 
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problem  ,  the  results  for  the  shock  trajectory  are  almost  identi¬ 
cal  ^  It  should  bn  noted  that  in  the  ncn-similar  so- 

(17) 

lution  of  Oshima  ,  the  mass  integral  is  not  even  satisfied. 
However,  by  requiring  that  the  total  energy  is  conserved,  Oshima' s 
solution  for  the  shock  trajectory  is  quite  accurate  when  compared 
to  the  more  rigorous  solution  of  Sakurai'  whore  both  the  mass 

and  energy  are  conserved. 

Based  on  the  insensitivity  of  the  shock  trajectory  on  the 
flow  distributions,  we  therefore  ignore  the  influence  of  the  chemical 
reactions  on  the  hydro  , mmic  flow  structure  and  assume  the  density, 
particle  veloc  ity  .  •'  ne  pressure  profiles  behind  the  spherical  wave 


lift,*)  *  *s(&) 


where  f!,~  ,  denote  the  first  and  second  derivative  of  with 

respect  to  time.  The  above  forms  for  the  flow  distributions  are 

(  1 

essentially  those  for  a  non-reacting  blast  wave  '  and  they  satisfy 


for  the  density 


the  conservation  integrals.  The  exponent 
profile  given  by  Eq .  2.3  can  be  evaluated  from  the  mass  integral 
(i.e.,  Eq .  2.1)  as 


f(4)  - 


(2.6) 


where  the  subscript  "  j  "  denotes  the  condition  at  the  shock  front. 

By  assuming  the  flow  structure  to  be  as  specified  by  Eqs.  2.3  to  2.5, 
we  no  longer  have  to  concern  ourselves  with  the  hydrodynamic  conserva¬ 
tion  equations. 

With  the  induction  zone  thickness  d(4)  finite,  the  total 
chemical  energy  released  in  the  flow  at  any  instant  will  be  given  by 


the  integral 


The  second  approximation  used  in  the  present  model  is  to  consider  the 
shock  as  reactive,  but  with  chemical  energy  Qe(£)  per  unit  mass 
released  instantaneously  at  the  front  instead  of  the 

total  chemical  energy  per  unit  mass  of  the  mixture  (j D  in  the  J-T-Z 
and  R-B-W  models.  This  effective  chemical  energy  Q,&)  depends 
on  the  local  induction  zone  thickness  and  the  dependence  can  be  obtained 
by  considering  the  conservation  of  chemical  energy  at  any  instant  of 


time,  i.e., 


4rcr}'gQe0t)dr 


Using  Eqs.  2.3  and  2.6,  the  effective  chemical  energy  Q<?  J 


can 
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bo  obtained  as 


QM)  -- 


/  -  c—'  i 

1  m) ' 


(2.7) 


If  the  induction  zone  thickness  d Cl )  Q  t  then  ^ 

and  the  present  model  reduces  to  the  R-B-W  model  i or  a  discontinuity 
detonation  front.  The  coupling  between  the  sb'ck  motion  and  chemical 
reactions  is  essentially  described  by  the  re!  tionship  llq .  2.7.  The 
degree  of  coupling  depends  on  the  induction  zone  thickness  dm . 

With  the  effective  chemical  energy  QCM)  released  at 
the  shock  front  instantaneously,  t he  boundary  conditions  can  be  written 
as 


/? 


W 


Y-s 

1  / 


(2.8) 


U,(£U)  = 


/-/■s-y 


y-f-i 


(2.9) 


(?,(%,■£)  = 

/?  R>) 


where 


y +ys> 

reno 


s  -  lO-yf-k yf 


'L 


(2.10) 


(2.H) 
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li--. 


(2.12) 


=  A  =  «*■  = 


7 


y?  4s 


(2.13) 


Note  that  there  are  two  roots  for  from  Eq.  2.11  with  the  positive 
root  representing  the  overdriven  detonation  solution  and  the  negative 
root  for  the  weak  detonation  solution.  The  unique  solution  for  ^  =  Q 
corresponds  to  the  C-J  solution.  The  positive  roe"  for  £  is  used  in 
the  present  analysis  since  we  are  concerned  with  the  decay  of  an  ini¬ 
tially  overdriven  wave  to  its  final  steady  state  velocity  when  S  -  O  • 
Since  Clc(x J  is  a  function  of  the  shock  radius,  the  final  steady 
state  velocity  depends  on  the  solution  and  the  magnitude  of 

QeC^)  is  determined  when  ^  -  O  . 

The  third  approximation  in  the  formulation  of  the  present 
model  is  the  specification  of  the  functional  form  for  the  induction 
scone  thickness  m  .  From  the  results  of  the  previous  study  of 
the  finite  kinetic  rate  blast  model  ^  and  experimental  observa¬ 


tions,  the  following  form  is  chosen 


dU)  *  A(l~e 


-X2' 


(2.U.) 


where 


yc  =  Q. 


(2.15) 


Sr‘  Rs/fZo 


(2.16) 


and  Ro  is  the  characteristic  explosion  length  commonly  used  in 
blast  wave  scaling  defined  as 


o  .  /£_ 

t-^/g 


(2.1V) 


3 

Writing  the  exponent:  'X-i  in  Eq.  2.14  in  terms  of  the  fundamental 
parameters,  we  see  that 


3 

'Xi  " 


%7TJ°QR? 


which  is  simply  the  ratio  of  the  total  chemical  energy  to  the 

3 

initiation  at  any  instant  of  time.  When  9^2-  /  ,  the 

# 

shock  motion  is  dominated  by  the  initiation  energy  and  the  front 
is  essentially  a  strong  shock  and  cfC£ )  — ^  O  •  When  2.  / 

the  chemical  energy  dominates  the  shock  motion  and  oY-/ )  — >-  /] 
which  is  some  constant  depending  on  the  particular  magnitude  of  final 
steady  state  shock  ‘"trength  of  the  solution,  both  these  limiting  be¬ 
haviours  as  well  as  the  exponential  nature  of  <d("t)  given  by  Eq. 

2.14  correspond  to  experimental  observations. 

As  a  final  step  in  the.  formulation  of  the  model,  a  lower 
limiting  shock  strength  where  chemical,  reactions  are  completely  de¬ 
coupled  from  the  shock  motion  has  to  be  specified.  This  is  based  on 
the  existence  of  a  temperature  limit  below  which  auto- ignit ion  of  the 
gas  behind  the  shock  is  no  longer  possible.  This  limit  is  ill- 
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defined  experimentally,  but  for  most  practical  gaseous  explosives 

o' V" 

of  interest,  the  limiting  temperature  is  of  the  order  of  1000  K. 

In  the  present  study  we  shall  take  a  critical  shock  strength 
for  the  lower  limit.  In  other  words,  if  the  shock  decays  below, 
in  any  particular  solution,  we  will  completely  decouple  any  chemical 

V 

reactions  from  the  shock  motion.  Hence  we  set  the  limiting  condit'iops 


for  Eq.  2.7  as 


Qe  =  QFfy*) 


mo 


.(2.18) 


=  O 


It  should  be  noted  that  the  "cut-off"  specified  in  Eq .  2.18  is  an 
abrupt  step-function  type  of  behaviour  while  in  reality  it  is  a 
smooth  function.  However,  due  to  the  sharp  exponential  increase  of 
dtt)  near  the  auto- ignition  limit,  the  assumed  step  function 
behaviour  is  quite  justified.  The  essential  steps  in  the  formulation 
of  the  present  model  have  now  been  specified. 


3.  Analysis 


For  convenience  in  the  analysis  and  numerical  computations, 
we  shall  recast  the  basic  equations  in  the  previous  section  in  the 
familiar  blast  wave  parameters  Replacing  the  independent 
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/ariables  (  r  ,t  )  by  (  7  )  where 

«wi) 

*,=  -  ft-  .  ^.1 


(3.1) 


r* 


/?.  * 


the  dependent  variables  (  >  y°  ;  &  )  by  (  -f  ,  ,  <£  ) 


where 


^  ■  H" 

•  4«' 

fo 


(3.2) 


i 

& 


The  density,  particle  velocity  and  the  pressure  profiles  behind  the 


spherical  wave  as  given  previously  by  Eqs.  2.3  to  2.3  become 


f  * 


(3.3) 


f 


(3.  A) 


*1  [  faj  •  '  y 


(3.3) 
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where 


Q-.  &&/£ 


(3.6) 


The  parameter  0  denotes  whether  the  shock  is  decelerating 
or  accelerating  &>o.  The  exponent  tci)  of  the  density 
profile  in  Eq.  3.3  can  be  written  in  terms  of  the  blast  wave  parameters 


as 


(3.7) 


The  subscript  /  in  Eqs.  3.3  to  3.5  denotes  the  conditions  immediately 
behind  the  shock  front  and  are  given  previously  by  Eqs. 

2.8  to  2.10. 

Substituting  the  profiles  for  f.f  and  j  into  the  energy 
integral,  we  obtain  the  following  equation 

t 

]L3Cr->) 

where  7C  ,  and  Ro  are  defined  previously  by  Eqs.  2,15  to  2.17 
respectively.  From  the  definition  of  0-  given  by  Eq .  3.6,  we  obtain 


(3.8) 


the  relationship 
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dy  = 

c/2  2- 


(3.9) 


so  that  the  term 


in  Eq.  3.8  may  be  written  as 


djjfit  _  c&i  +  -  c/d,  -  Jt-  2d, 

dr)  dy  £y  <2Qy 


(3.10) 


Evaluating 
at  the  front  for 


and  2>d,/z2- 

from  the  boundary  condition 
^  (i.e.,  Eq.  2.9),  we  obtain 


=  <gy3)#  ^  (t£> 

D#fiO  s  f.  Aft-%) 

and 


(3.11) 


(3.12) 


x  *,3  (3*13) 

=  (finyzsk .  (t-DiSz^Je  ) 

£*(/-%)  s+ 

^  ar.d  ^  are  defined  previously  by  Eqs.  2.11  and  2.12  respectively 
wh<Vc  and  ^  are  defined  as 
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s--  d(i)/& 


(3.14) 


Using  Eq.  3.10,  we  can  solve  for  Q  from  Eq.  3.8  as 


e  -  x&xpaz  )l0*  d'-if) 


-A  +■  yt<tCf-<$)  _ 

2(7-0  W3) 


(3.15) 


iCst/gtDs] 

Substituting  the  above  equation  into  Eq.  3.8,  we  obtain  a  differential 
equation  for  the  change  of  the  shock  strength  *j  with  the  radius  2?  , 


i.e., 


P  '■ 

d? 


(3.16) 


which  may  be  integrated  numerically  using,  for  example,  the  standard 
Runge-Kutta  method. 

To  start  the  numerical  integration  wo  must  obtain  the  value 
of  y at  some  small  initial  value  of  .  From  Eq .  3.8,  we  note  that 
as  £-5*0  ,  y  ->  2 3  since  the  left  hand  side  of  the  equation 
is  finite.  Hence  the  form  of  the  expansion 


for  small  values 


Hence  the  shock  behaves  like  a  strong  point  blast  wave. 
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For  very  large  initiation  energy,  00  or 

•fig,  — ^  and  hence  from  Eqs.  2.14  and  2.7,  dte)-*»o  so  that 

— *■-  •  1°  other  words,  the  full  chemical  energy  per  unit 

mass  is  released  instantaneously  at  the  front  and  the  present  solution 
recovers  the  reacting  blast  wave  model  Alternatively,  we 

can  assume  the  front  to  be  discontinuous  by  simply  taking  /^  =  Q 
and  again  the  reacting  blast  wave  model  is  recovered.  However,  in  this 
case,  the  initiation  energy  can  be  arbitrary  and  in  the  limit 
we  recover  the  J-T-Z  model. 

Taking  C£=  O  ,  the  boundary  conditions  (i.e.,  Eqs.  2.8  to 
2.10)  reduce  to  the  ordinary  Rankine-Hvgoniot  equations  for  a  normal 
shock  in  a  perfect  gas.  The  present  solution  then  becomes  the  solution 

/  J  o\ 

for  a  non-reacting  blast  wave  .  With  the  non-reacting  hydrodynamic 
flow  structure  determined,  the  chemistry  of  the  reactions  between  the 
shock  can  be  computed.  This  in  essence  is  the  finite  kinetic  rate 
model  studied  by  the  Authors  and  Lundstrom  and  Oppenheim  , 


Therefore  we  see  that  the  present  model  is  consistent  with  and  unifies 
all  the  previous  theoretical  models. 
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4.  Results  and  Discussions 


For  an  explosive  mixture  at  given  initial  conditions  we 
shall  specify  the  following  parameters:  the  specific  heat  ratio  y  , 
the  initial  pressure  and  sound  speed  C0  ,  the  planar  C-J  Mach  No. 
Mct  (hence  the  chemical  energy  per  unit  mass  Q),  and  the  critical 
shock  Mach  No.  Mc  corresponding  to  the  auto-ignition  limit  of  the 
mixture.  For  the  system  a  A  -  4  studied  in  this  paper,  we 
choose =  £  which  corresponds  to  an  auto- ignition  temperature  limit 

Tc  ioooV 

.  For  any  chosen  value  of  the  initiation  energy  parameter  <c  ~ 

,  we  can  obtain  the  solutions  for  the  shock  trajectory  ■m) 
and  the  variation  of  shock  strength  with  radius  J(2)  by  integrating 
Eqs.  3.16  and  3.20.  To  start  the  numerical  integration,  the  perturbation 
expressions  given  by  Eqs.  3.17  and  3.18  are  used  to  evaluate  the  starting 
values  for  ^  and  Q  for  some  small  chosen  value  of =  z?o  .  »<>  the 
present  work  2*>  was  chosen  to  be  0,05  and  the  integration  was  found  to 
proceed  smoothly  for  increasing  values  of  £.  without  oscillation.-.. 

It  is  important  to  consider  first  the  general  properties  of 
the  solutions  for  )  ,  ^(-i)  and  SCi)  for  various  values 

of  the  initiation  energy  parameter  G  .  A  typical  example  for  the  case 
of  stoichiometric  CzfJi  -  Oj.  undiluted  mixture  at  an  initial  pressure 
^9  =  100  torr  and  temperature  =  300  K  is  shown  in  Fig.  1.  Tne 

solutions  for  &(■£ )  ■  jr  various  values  of  €  all  start  out  at  the 

limiting  strong  shock  value  of  @Co)  =  -3/2  and  increase  towards  zero 
as  the  shock  expands  and  the  influence  of  the  chemical  energy  release 
becomes  progressively  important  with  larger  bhock  radius  .  The 
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solutions  of  s&)  all  start  out  at  its  initial  value  SCo)  =  /  , 
which  corresponds  to  the  strong  shock  limit  9?  =  O  (i.e.,  ~  ) 

and  decrease  monotomically  to  zero  as  the  shock  expands  and  decelerates. 
From  Kq .  3.15,  we  see  that  when  S  -  O  >  the  parameter  &  =  O  ,  which 
implies  that  the  shock  acceleration  dRjcH  ,  /$  --  O  .  Hence  as 
S O  ,  the  shock  ceases  to  decelerate  and  will  propagate  there¬ 
after  at  a  final  steady  state  ve.locity  fiff  ,  Since  corresponds 

to  S  -  &  >  we  sec  rrom  the  boundary  conditions  (i.e.,  Eqs.  2.8  to 
2.10)  that  this  solution  is  one  where  the  Payleigh  line  is  tangent  to 
the  equilibrium  Hugoniot  curve  for  the  particular  value  of  Qe  obtained . 
Hence  this  final  steady  state  velocity  can  be  interpreted  as  the 

local  C-J  velocity.  As  Mj*"  is  a  function  of  <r  ,  a  different  final 
steady  state  velocity  is  obtained  for  different  values  of  the  initiation 
energy  parameter  £  .  This  implies  that  for  spherical  detonations,  the 
final  steady  state  C-J  velocity  depends  on  the  magnitude  of  the  initiation 
energy  E0  .  The  present .conclusion  is  then  contradictory  to  the  classi¬ 
cal  theory  where,  for  a  given  explosive,  the  steady  state  C-J  velocity 
depends  only  on  the  properties  of  the  explosive  itself. 

From  Fig.  1,  it  is  seen  that  for  small  values  of  £  (i.e., 
large  initiation  energy)  the  solutions  for  0(2 )  and  S(2)  indicate 
that  ec*)  o  simultaneously  as  .  However, 

for  larger  values  of  £  ,  we  note  that  &(•£)  -*’•  G>  when  Sfe) 
is  still  finite  and  9(*)  hence  becomes  positive  and  continues  to 
increase  as  Sfr)  monotomically  decreases  to  zero.  When  SCi)  =  o, 
the  value  of  Q( 2)  becomes  multi-value  since  from  Eq .  3.15,  @  must 

vanish  when  S  -  Q  .  Physically  this  implies  that  the  shock  over-deco- 
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lerates  at  first  to  a  minimum  velocity  and  as  S&)  is  still 
finite,  the  shock  accelerates  again  until  SOz)  =  o  when  the 
acceleration  must  suddenly  cease  (i.e.,  0  ~  O  )  ^or  t^ie  shock  to 

propagate  subsequently  at  a  steady  state  velocity  •  The  dis¬ 

continuous  jump  in  the  shock  acceleration  is  illustrated  in  the  large 
abrupt  change  in  the  solution  for  Tor  the  case  of  G  =  0.13 

shown  in  Fig.  1.  This  sudden  transition  demonstrates  a  type  of  in¬ 
stability  which  occurs  when  the  hydrodynamic  flow  structure  is  not 
coupled  smoothly  to  the  chemical  reactions.  It  should  be  noted  that 
in  the  numerical  studies  of  the  deceleration  of  an  initially  planar 
overdriven  detonation  wave  ^  ^  the  onset  of  instability  is  also 
observed  as  the  shock  approaches  its  final  steady  state  C-J  condition. 
This  type  of  instability  is  inherent  in  all  detonation  problems  and 
within  the  context  of  the  present  model,  more  detailed  information 
on  this  instability  mechanism  cannot  be  obtained.  However,  for  small 
values  of  the  initiation  energy  parameter  ,  the  ampl  .aides  of  the 
jumps  in  Bti)  are  sufficiently  sma  1  and  are  within  the  round  off 
errors  of  the  finite  step  sizes  of  the  numerical  integration  scheme. 
Hence  we  shall  assume  that  for  small  values  of  (i.e.,  O  £  €  &  ..05 
for  the  particular  case  shown  in  Fig.  1)  the  deceleration  of  the  sphe¬ 
rical  detonation  wave  to  its  final  steady  state  velocity  is 

smooth . 


The  stable  behaviour  of  the  solutions  for  small,  values  of  (r 
can  best  be  illustrated  by  the  solutions  for  .  We  note  that 

Cor  the  range  of  O /  <  o'. .  the  shock  strength  decays  smooth¬ 

ly  from  its  initial  value  vj(o)  -  Q  to  some*  final  steady  state  value 
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-  .  Outside  this  range  05€  ?£  .05,  the  shock 

decays  to  a  minimum  and  accelerates  again.  The  transition  to  its 

final  steady  state  value  ^^"when  5=0  is  discontinuous.  This 

is  illustrated  clearly  for  the  case  £  =  0.13.  Since  for  each  shock 

radius  it  the  local  value  of  6)^  can  be  evaluated  from  Eqs.  2.7  and 

2.14,  we  can  determine  the  local  C-J  velocity  by  taking  <£  =  C>  and 

solving  for  y\  from  Eq.  2.11.  The  variations  of  this  local  C-J  value 

A,  ,  / 

of  'rj  ( t-y  with  radius  for  the  various  values  of  £  are  also  shown 
in  Fig.  1.  We  note  that  the  corresponding  yf-2i)  is  always  less  than 
the  local  C-J  value  of  y'te)  during  the  decay  (the  shock  speed  is 
always  greater  than  its  local  C-J  velocity).  Hence  the  detonation 
wave  is  overdriven  at  all  times  until  5=0  when  r>jC?-)  =  = 

j  and  the  final  steady  state  velocity  is -obtained. 

In  Fig.  2,  the  variation  of  the  detonation  Mach  No. 
with  radius  <£  =  ^  //?o  is  shown.  For- small  shock  radii  (i.e,, 

)  where  the  initiation  energy  dominates  the  shock  motion, 

all  the  curves  coalesce  and  the  decay  is  essentially  that  of  a  strong 

,  IL 

blast  (i.e.,  f»s  /  £  ).  For  example,  all  the  solutions  for 

are  indiscernible  from  each  other  for  ^  0.2,  which 
corresponds  to  the  energy  ratio 


\S  Cr  „ oe 


£q 


In  other  words,  the  total  chemical  energy  enclosed  by  the  wave  is 
about  8$  of  the  initiation  energy  (E0  on 1 y .  Significant  departures 

t 

from  the  blast  wave  solution  should  occur  when  the  energy  ratio 
which,  for  the  particular  ease  shown,  is  at  a  shock  radius  2  ^  .46 
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where  the  differences  between  the  various  solutions  are  quite  pro¬ 
nounced.  The  final  steady  state  conditions  (i.e.,  Mi  /%»*  > 

correspond  to  a  shock  radius  ^  Cr*  0.6.  This  yields  an  energy 

3 

ratio  a\.Z.  •  In  other  words,  when  the  chemical  energy  fi¬ 

nally  dominates  the  shock  motion, 

4-rc^CiRs 


Although  the  so-called  smooth  stable  regime  Q~<~  0.5 

for  the  decay  to  the  final  steady  state  velocity  is  arbitrarily  set 
based  on  the  vanishingly  small  abrupt  change  of  slope  as  /V?g\ 

it  is  interesting  to  note  that  this  stable  limit  is  in  accord  with  the 

(19) 

stability  criterion  of  Shchelkin  .  Using  the  final  steady  state 

velocity  =5.7  for  the  limiting  value  of  (~  -  0.5,  the  induction 

time  computed  using  White's  data  for  this  mixture  (i.e., 

J  V  -/&£/•?■  /fJco/i/SB'D  is  r^=  .056  x  10“6 

_  -5 

sec.,  wDich  is  about  twice  the  value  LCj-  =  .029  x  10  sec.  for  the 
planar  C-J  velocity  of  =  7.015.  According  to  Shchelkin,  if  a 

random  disturbance  results  in  an  increase  in  the  induction  time  by  an 
order  of  magnitude  equal  to  or  greater  than  the  initial  delay  time  it¬ 
self,  then  the  detonation  front  loses  its  stability.  For  spherical 
detonations,  the  increase  in  the  induction  time  is  due  to  the  lateral 
expansion  associated  with  the  three-dimensional  geometry. 

For  <£-  ^  .05  the  detonation  wave  decays  to  a  minimum 

and  then  accelerates  again.  The  transition  to  its  final  steady  state 
vo loc ity  when  3  C>  is  discontinuous.  This  corresponds  to  the 


exper iisoncal  observations  of  the  propagation  of  a  spherical  detonation 
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wave  in  Che  critical  initiation  energy  regime  ,  Tlie  detonation 
decays  continuously  to  a  velocity  corresponding  to  the  auto-ignition 
limit  of  fll(2,  2*  3  .  The  shock-reaction  zone  complex  then  pro¬ 
pagates  at  this  minimum  velocity  for  a  short  period  of  the  order  of: 

10  pscc.  Tiie  onset  of  instability  is  observed  in  the  form  of  the 
appearance  of  local  ignition  centers.  The  rapid  development  of  these 
local  hot  spots  leads  to  a  sudden  transition  to  a  highly  asymmetrical 
multi-headed  wave.  From  the  experimental  record  of  our  previous  pa¬ 
per  which  was  for  stoichiometric  02~^2**2  at  torr,  the  ini¬ 
tiation  energy  =  0.3  joules  and  the  value  of  /\  ~  2  mm  yielding 

a  value  of  .18.  From  Fig.  2  we  note  that  this  corresponds 

closely  to  the  critical  value  of  -  0.13  for  direct  initiation  to 
be  possible. 

For  very  small  initiation  energies  (i.e,,  £  -  ^/f?0  is 
large),  the  wave  decays  ,  be  low  the  specified  limiting  strength  for 
auto- ignition  (i.e.,  f£c  =  3).  The  shock  completely  decouples  from 
further  chemical  reactions  and  decays  asymptotically  to  an  acoustic 

wave/!'/-  =  1.  This  regime  is  the  so-called  sub-critical  energy 

,  ..  ,  .  .  (10) 

regime  defined  in  a  previous  paper  . 

From  the  present  theory  the  critical  energy  for  direct 
initiation  can  be  determined  since  the  critical  value  of  Lhe 
initiation  energy  parameter  is  known.  To  eva¬ 

luate  ,  we  first  compute  4,#/  using  the  formula  of  White 

-Z0.8/  *-!?3eo/4'Ser  ) 

based  on  the  final  steady  state  velocity/^  *■  3.  With  ) 

known,  the  critical  energy  can  be  determined  as  follows: 
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&c  f  is-  f 


hence 


Ec  "(Z^y^ 

The  dependence  of  £^,on  as  well  as  Up.)  from  White's 
formula  for  the  two  cases  of  stoichiometric  and  equi-molar  C2H2~^2 
mixtures  arc  illustrated  in  Fig. 3.  The  experimental  points  for  the 

i 

narrow  pressure  range  .04  ^/4  ^  .15  atm.  are  taken  from  our 

(10) 

previous  paper  where  a  laser  spark  of  20  Msec,  duration  is  used 

as  the  energy  source.  The  agreement  is  quite  reasonable  in  view  of 
the  uncertainty  in  the  use  of  the  kinetic  data  near  the  low  tempera- 
ture  limit.  It  should  be  noted  that  although  more  data  are  available 
for  the  critical  energy  for  direct  initiation  for  these  explosive 
mixtures,  they  cannot  be  used  to  compare  with  the  present  theoretical 
results.  It  was  pointed  out  in  a  recent  paper  that  the  critical 

energy  can  differ  by  a  few  orders  of  magnitude  depending  on  the  total 
discharge  time.  In  order  to  compare  with  the  present  theoretical 


Be 


00  ),  these 


results  based  on  a  point  blast  model  (i.e., 

&{ 

existing  data  must  be  reduced  to  the  limiting  value  corresponding  co 
an  infinite  power  density  of  the  ignition  source.  Unfortunately,  the 
details  of  the  discharge  characteristics  for  these  experiments  where 
the  data  are  reported  have  not  been  stated  to  permit  the  appropriate 
data  reduction.  The  laser  spark  is  perhaps  the  most  powerful  point 
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energy  source  available  which  approximates  closest  to  the  theoretical 
requirement  of  point  blast  theory. 

To  confirm  the  existence  of  a  final  steady  state  velocity 

for  spherical  detonations  that  are  below  the  classical  value  for  the 

explosive,  a  comparison  between  the  present  results  and  the  recent 

(21)  (22) 

experimental  data  of  Brossard  and  Struck  is  made.  In  Bros- 

sard's  experiments,  a  1  meter  diameter  spherical  bomb  is  used  and  the 
detonation  velocity  is  measured  by  the  microwave  doppler  technique  to 
an  accuracy  of  0.4$.  Various  mixtures  (C^H^,  (^H^,  ^2^2  °2^ 

at  an  initial  temperature  of  293°K  and  .1  ^ ^  2.8  bar  have  been 

studied.  In  view  of  the  large  diameter  chamber  used  to  ensure  steady 
state  conditions  as  well  as  the  accuracy  of  the  wave  velocity  measure¬ 
ment  by  the  doppler  shift  method,  Brossard 's  results  represent  the 
most  reliable  existing  data  on  spherical  detonation  velocities.  In 
S truck's  experiments,  a  much  smaller  chamber  of  100. mm  diameter  was 
used.  The  detonation  wave  velocity  was  measured  using  high  speed  mo¬ 
vie  schlieren.  The  present  theoretical  results  for  e‘h  > 

are  shown  in  Fig.  4  for  both  stoichiometric  and  equi-molar 
mixtures.  In  order  to  compare  with  Broward  and  Struck'c  experimental 
data,  the  variations  of  with  initial  pressure  ^  ^or  constant 

initiation  energies  must  be  obtained.  To  derive  these  constant 
energy  velocity  curves,  we  proceed  as  follows:  Suppose  we  desire  i 
for  a  particular  energy  .  We  note  first  of  all  that  this  curve 
must  originate  at  the  critical  condition  where  e0-  ec.  From  Fig. 
3,  the  value  of  and  can  be  obtained  for  any  value  of  • 
Hence  at  any  pressure  j?  ,  wo  may  write 
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AcCf>) 1 A  } 


For  stoichiometric  C^-O,,,  £c  *  0.15  while  for  equi-molar  mixtures, 
^C.  =  0.13.  Since  for  various  values  of  £  are  determined  from 

the  solution,  )  of  any  particular  value  of  fi^can  be  derived. 

From  fig.  4  we  note  that  the  agreement  between  the  present  theory  with 
the  data  of  Brossard  and  Struck  is  quite  good.  It  should  be  pointed 
out  that  from  the  present  theory  the  detonation  velocity  depends  only 


weakly  on  the  initiation  energy  in  the  high  pressure  range  .55  £ j  ^  3. 
Therefore  although  the  initiation  energy  of  both  Brossard  and  Struck 
has  not  been  properly  reduced  to  the  limiting  value  corresponding  to 
infinite  power  density  as  required  by  the  point  blast  model,  the  agree¬ 
ment  is  still  very  good. 

Theoretical  attempts  ^  2  2^  have  been  made  to  predict 

the  cellular  dimension  of  a  multiheaded  detonation  front.  Using  the 
present  theory,  the  longitudinal  dimension  of  a  "detonation  cell"  can 
be  evaluated  if  it  is  assumed  that  the  motion  of  a  "detonation  wavelet" 
in  between  collisions  of  transverse  waves  behaves  like  a  decaying  react¬ 
ing  blast  as  described  by  the  present  theory.  Based  on  the  experimental 
evidence  that  the  minimum  strength  of  a  "detonation  wavelet"  just  prior 
to  the  re -energization  by  the  collision  of  a  pair  of  transverse  waves 
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is  about  /p$  ^  3  which  corresponds  to  the  critical  limiting 
strength  for  auto- ignition,  we  shall  take  the  shock  radius  •£.c  for 
the  case  where  £  =  <£c  as  the  longitudinal  dimension  of  the  detona¬ 
tion  cell.  From  F '$•  2,  we  note  that  £c  Or  .6  and  the  value  of 

S*  .15.  These  values  (i.e.,  ^  )  are  almost  the  same 

for  a  wide  pressure  range  for  stoichiometric  mixtures  while 

for  equi-molar  mixtures  (£c  >  .13  and  ^  ^  .6.  Using  these 

values,  the  longitudinal  cell  dimension  Fj  can  be  obtained  through 
the  relationship 

n  -  ic  AcCfe) 
d  e 

Using  White's  formula  to  determine  /S.cCjl)  >  the  variation  of  Rj 

with  initial  pressures  hQ  can  be  obtained.  To  compare  with  the 

(2b\ 

experimental  data  of  Strehlow  '  ,  the  longitudinal  dimension  must 

be  multiplied  by  a  factor  of  about  0.62  based  on  the  definition  he 
used  for  the  transverse  wave  spacing.  This  comparison  is  shown  in 
Fig.  5  and,  as  can  be  observed,  the  agreement  is  reasonable.  It  should 
be  pointed  out  that  there  is  no  concrete  experimental  evidence  to  sup¬ 
port  the  assumption  that  the  motion  of  a  detonation  wavelet  can  be  mo¬ 
delled  by  a  point  blast  wave;  In  fact,  as  Strehlow  pointed  out,  the 
center  of  curvature  of  the  detonation  wavelet  does  not  corespond  to  the 
point  of  collision  of  the  transverse  wave  where  the  "decaying  wavelet" 
first  originates.  Also,  the  impulsive  release  of  energy  by  the  collision 
of  transverse  waves  is  by  no  means  instantaneous,  lienee  the  subsequent 
decay  of  the  "wavelet"  does  not  correspond  to  that  from  the  present  model. 
However,  in  view  of  the  similarity  in  the  gross  features  of  the  decay 
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processes  between  experiments  and  the  present  theory,  it  seems 
reasonable  to  adopt  the  blast  wave  model  as  a  first  approximation. 

5 .  Conclusions 

It  should  be  emphasized  that  the  most  important  results 
of  the  present  theory  are  that  it  predicts  the  existence  of  steady 
state  sub-Chapman-Jouguet  velocities  which  depend  on  the  details  of 
the  initiation  processes.  In  other  words,  the  final  steady  motion 
of  the  detonation  front  bears  the  memory  of  the  initiation  condition 
through  the  hydrodynamic  flow  structure  generated  by  the  ignition 
source.  This  is  supported  by  the  recent  accurate  experimental  mea¬ 
surements  made  by  gross a rd .  If  the  present  theory,  together  with 
Brossard's  as.  well  as  Struck's  experimental  results,  are  to  be  taken 
as  correct,  then  the  classical  Chapman-Jouguot  theory  of  detonation 
waves  must  be  modified.  The  abundance  of  the  experimental  results 
gathered  in  the  past  decade  strongly  indicates  the  inadequacy  of  the 

Chapman-Jouguet  theory.  The  recent  verbal  critique  of  the  Chapman- 

( 

Jouguet  theory  by  Professor  Hanson  at  the  13th  Combustion  Symposium 
emphasizes  the  importance  of  the  influence  of  the  hydrodynamic  flow 
structure  in  the  wake  of  the  expanding  detonation  on  its  motion.  The 
present  theory  can  be  considered  as  a  particular  case  where  the  hydro- 
dynamic  flow  structure  is  peculiar  to  that  of  a  point  energy  source. 

It  would  be  important  to  extend  the  present  model  to  other  types  of 
hydrodynamic  flow  structure  in  order  to  establish  whether  the  predic¬ 
tions  from  t he  present  model  are  indeed  universal.  Equally  important 
will  be  more  careful  experiments  on  spherical  detonations  in  the  low 
pressure,  low  initiation  regime  where  the  effects  of  the  experiments 
can  be  manifested  in  a  more  pronounced  fashion. 
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Captions 

The  solutions  of  ,  9(i)  >  S2(2)  and  %  (*) 

for  various  values  of  6  for  stoichiometric  in*xtures 

at  "]>0  =  100  torr,  "To  =  300  K. 

The  variation  of  shock  strength  with  shock  radius  for 
various  values  of  6  (stoichiometric  C2H2-02,  =  100  torr, 

T0  =  300°K). 

The  variation  of  the  critical  energy  Ec  for  direct  initiation 
and  the  critical  induction  zone  thickness  Ac  with  initial 
pressures  for  stoichiometric  and  equi-molar  02H2~02  mixtures, 
Comparison  of  the  detonation  velocities  in  stoichiometric  and 
equi-molar  C2H2-02  mixtures  with  the  experimental  data  of 
grossard  (Ref.  21)  and  Struck  (Ref,  22). 

Comparison  of  the  predicted  transverse  wave  spacings  with 
the  experimental  data  of  Strehlow  (Ref.  26), 


02-C2H2  mixtures 


